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Abstract. We have performed 3D radiation MHD simulations extending from the convec- 
tion zone to the corona covering a box 16 Mm 3 at 32 km spatial resolution. The simulations 
show very fine structure in the chromosphere with acoustic shocks interacting with the mag- 
netic field. Magnetic flux concentrations have a temperature lower than the surroundings in 
the photosphere but higher in the low chromosphere. The heating is there mostly through 
ohmic dissipation preferentially at the edges of the flux concentrations. The magnetic field 
is often wound up around the flux concentrations. When acoustic waves travel up along the 
field this topology leads to swirling motions seen in chromospheric diagnostic lines such as 
the calcium infrared triplet. 
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1. Introduction 

Inferring physical conditions in the solar chro- 
mosphere is difficult, because the few spectral 
lines that carry information from this region 
(e.g. hydrogen Balmer-a, Call H and K lines 
and infrared triplet and He 1083 nm line) are 
formed outside local thermodynamic equilib- 
rium (LTE). Thus they do not directly couple 
to the local conditions. In addition, the chro- 
mosphere is very dynamic with large varia- 
tions in temperature and density and a concept 
like "formation height" is not very useful. We 
also go from plasma domination in the pho- 
tosphere to a magnetically dominated regime 
in the upper chromosphere. This means that 



the chromosphere can sustain a multitude of 
wavemodes that become degenerate and cou- 
ple when the sound speed equals the Alfven 
speed. Major progress in our understanding is 
dependent on comparison between detailed ob- 
servations and detailed modelling. Such mod- 
elling needs to go all the way from the convec- 
tion zone (where waves are excited and mag- 
netic field foot-point braiding takes place) to 
the corona (to enable magnetic field connec- 
tivity) and include a reasonable treatment of 
the radiative tra nsfer. The first su ch modelling 
was reported by lHansteeril (120041) . These sim- 
ulations have been used to study flux emer- 
gence through the photosphere into the corona 
dMartmez-Svkora et al]l2008[l2009al) a nd driv- 
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ing of spicules (Martmez-Sykora et al. 2009b) 
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We here analyze a high resolution simulation 
focusing on the chromospheric structure, heat- 
ing and dynamics close to magnetic flux con- 
centrations. The layout of this contribution is 
as follows: in Section|2]we outline the methods 
used, in Section[3]we describe how the simula- 
tion analyzed was set up, in Section [4] we give 
the results and we finish in Section[5]with con- 
clusions. 

2. Methods 

Earlier simulations were carried out with the 
"Oslo Stagger Code", a code where the par- 
allelization was performed assuming shared 
memory and therefore restricted to rather few 
CPUs, see lHansteen. Carlsson. & Gud iksen 
(2007) for a description. We have now com- 
pleted a full rewrite using the Message Passing 
Interface (MPI) parallelization library and a 
completely new radiation transfer module em- 
ployi ng domain decomposition (IHavek et al.1 
1201 Oh . This new code, named BIFROST, is 
outline d below. A det ailed writeup is in prepa- 
ration dGudiksen et alb . 

We solve the equations of radiation- 
magnetohydrodynamics using a compact, 
high-order, finite difference scheme on a 
staggered mesh (6th order derivatives, 5th 
order interpolation). The equations are stepped 
forward in time using the explicit 3 rd order 
predic tor-corrector procedure of Hvman 
dl979l) . modified for variable time steps. 
High-order artificial diffusion is added both 
in the forms of a viscosity and in the form 
of a magnetic diffusivity. The radiative trans- 
fer equation is solved along 24 rays using 
a domain-decomposed short-characteristic 
method. The opacity from millions of spectral 
lines is inc luded using mu lti-group opacities in 
four bins (Nordlund 1982) modified to include 
the e ffects of coherent scattering dSkartlienl 
2000). Conduction along the magnetic field- 
lines is solved implicitly using a multigrid 
method. In addition to the radiative energy 
exchange treated with the detailed solution 
of the radiative transfer (including scattering 
effects) we include optically thin radiative 
losses in the the corona and radiative losses 
from strong lines of hydrogen and singly 



ionized calcium. These latter losses dominate 
the radiative lo s ses fr om the chromosphere 
dVernazza et alj Il98ll) . We include these 
radiative losses without assuming LTE by 
using an expression 

<D = e(T)e- T f(T) 

where O gives the radiative losses per electron 
per atom from a given element (hydrogen and 
singly ionized calcium, respectively), 6(T) is 
the sum of all collisional excitations (obtained 
from atomic data), T is temperature, r is an op- 
tical depth and e~ T /(r) is an effective escape 
probability. The optical depth r is set propor- 
tional to the vertical column mass with the pro- 
portionality constant obtained from the ID de- 
tailed non-LTE radiation hydrodynamic simu- 
lations of Carlsson & Stein (1992, 1995, 1997, 
2002) The function f(j) is obtained from a fit 
using the same ID simulations. 



3. Simulation setup 

The simulation analyzed in this contribution 
covers a region 16.6x16.6x15.5 Mm 3 with the 
bottom boundary in the convection zone 1.4 
Mm below the height where the mean optical 
depth at 500 nm is unity (in the following this 
height defines our z=0) and the top boundary 
in the corona 14.1 Mm above this height. Both 
the bottom and top boundaries are transparent. 
The simulation domain is periodic in both the 
horizontal directions with a fixed spacing of 
32.5 km while the vertical grid size is 24-28 
km up to 4 Mm height, increasing gradually 
above this height to a final maximum value of 
213 km at the top of the computational domain. 
The final simulation has 512x512x325 grid- 
points. The initial relaxation of the convection 
was performed at lower resolution. 

The mean, unsigned field strength at z=0 
is initially 40 G but falls to about 10 G at the 
beginning of the high-resolution run due to the 
fact that no field is injected into the computa- 
tional domain and the existing, tangled field, 
decays slowly. 
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4. Results 

After the increase of the resolution, there is a 
slow increase in the spatial power at high fre- 
quencies. After some 400 s the evolution starts 
to diverge from the continuation of the low res- 
olution run — granules fragment more easily 
and the velocity power at small spatial scales 
increases. The mean unsigned field strength 
also starts to increase, presumably due to local 
dynamo action made possible at the higher res- 
olution, and the mean unsigned field strength 
ends up at 15 G around t= 1500s and stays 
at this level until the end of the simulation at 
t=1670s. The vertical magnetic field strength 
at z=0, t=1500s is shown in Fig.Q] As can be 
expected from the low average unsigned field 
of 15 G there are large regions with only weak 
field. There are several regions with one polar- 
ity dominating but in general, both polarities 
co-exist. This is even more evident if the im- 
age scaling is adjusted to emphasize the weak 
field. The magnetic flux concentrations have a 
field strength close to 1 kG. 




X TMnnl 

Fig. 1. Vertical magnetic field strength at z=0, 
t= 1500 s. The colour scale has been clipped at 500 
G to emphasize weaker field. The maximum field 
strength is ± 1 kG. The square marks the region of 
interest shown in more detail in Fig. [3] 

The temperature structure of these mag- 
netic flux concentrations compared with the 



temperature structure of granules and inter- 
granular lanes is shown in Fig. [2] The average 
temperature for granules has been taken over 
all columns with a vertical upflow larger than 
100 ms 1 at z=0. The intergranular lane av- 
erage has been taken over all columns with a 
vertical downflow larger than 100 ms _I and 
the average for magnetic flux concentrations 
has been taken over all columns where the 
unsigned magnetic field in the photosphere is 
larger than 300 G. Intergranular lanes are hot- 
ter than granules above the surface (inverse 
granulation). Magnetic flux concentrations are 
cooler than intergranular lanes below the sur- 
face at a given geometric height but hotter than 
both granules and intergranular lanes in the up- 
per photosphere and lower chromosphere. 
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Fig. 2. Mean temperature as function of height for 
granules (dashed), intergranular lanes (dash-dotted) 
and magnetic flux concentrations (solid). 



This increased temperature in the upper 
photosphere-lower chromosphere of magnetic 
flux concentrations is clearly shown in Fig. [3] 
The region shown with a square in Fig. Q] 
is shown in more detail at a height of 0.29 
Mm with temperature, magnetic field strength, 
Joule heating and viscous heating displayed. It 
is clear that most magnetic flux concentrations 
have a higher temperature than the surround- 
ings at this height and that Joule heating takes 
place at the locations of high magnetic field 
strength. These locations evolve only slowly 
in time while the viscous heating is concen- 
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Fig. 3. Temperature (upper left panel), magnetic field strength (upper right), Joule heating (lower left) and 
viscous heating (lower right) at a height of 0.29 Mm at t= 1 500 s. At this height, magnetic flux concentrations 
are hotter than the surroundings due to Joule heating 



trated at the edges of granules evolving on 
much shorter timescales. 

Due to foot-point motions, the magnetic 
field topology in the chromosphere changes all 
the time. In the upper photosphere and lower 
chromosphere the plasma j3 (ratio between 
gas pressure and magnetic pressure) is larger 
than unity and the magnetic field is far from 
force-free. In the upper chromosphere, /3 
drops below unity and the field approaches 



a force-free state. This interplay between 
field-braiding and expansion into a force-free 
state often leads to spiral topologies around 
small magnetic flux concentrations. One 
example from the simulation is shown in Fig. [4] 
where the field-topology is shown around 
a magnetic flux concentration at t= 1010 s. 
Field-lines from the negative polarity in the 
flux concentration connect to various positive 
polarity patches with a clear spiral pattern as 
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Fig. 4. Magnetic field above a magnetic field concentration as seen from above (left panel) and in per- 
spective from the lower left corner (right panel). The vertical component of the field at z=0 is shown in 
greyscale. The sign of Bz is shown in colour with red field lines corresponding to negative Bz (dark in the 
plane at z=0). The highest field-lines shown extend to a height of 5.8 Mm. Note the spiral topology as seen 
from above and the abundance of low-lying almost horizontal field-lines. 
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Fig. 5. Synthesized narrow band image at the core of the chromospheric ,4854.2 nm line from singly ionized 
calcium at two times showing the apparent swirling motion. 



seen from above (left panel in Fig |4| and a 
large amount of almost horizontal (but curved) 
field (right panel in Fig |4|. At this instant in 
time in the simulation there is an acoustic wave 
that turn into a slow-mode wave in the low-/? 
region travelling along the curved field-lines. 



In synthesized, narrow band, images from the 
core of the /1854.2 nm line of singly ionized 
calcium (Fig|5]l this is visible as a swirling 
motion around the flux-concentration very 
similar to observations recently reported by 
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IWedemever-Bohm & Ro uppe van der Voor j 
d2009l) . 

5. Discussion and conclusions 

A multitude of phenomena are seen in these 
high resolution radiation-MHD simulations 
that extend from the convection zone to the 
corona and we have here only focused on 
the heating and dynamics around magnetic 
flux concentrations. Although these simu- 
lations are the most ambitious performed 
to date, there are still important ingredients 
missing. The current simulations do not 
include the effects of out-of-equilibrium 
hydrogen ionization even though it is known 
that such effects are very important for the 
energy balance of especial ly the middle- 
highe r chromosphere ([C arlsson & Steinl 
2002t Eeenaarts & Wedemever-Bohml 120061: 
Leenaarts et alj I2007I) . Including such effects 
would likely increase the amplitude of shocks 
in the middle chromosphere. Another limita- 
tion is the neglect of the back-radiation from 
the corona onto the upper chromosphere. This 
radiation would heat the chromosphere above 
about 1 Mm height dCarlsson & Steinl 12002b 
but have little effect on the heights discussed 
in this contribution. 
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